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ABSTRACT 

We present an analytical approach to the dynamical evolution of fireballs or axially symmetric jets expanding 
into an external medium, with application to Gamma-Ray Burst remnants. This method leads to numerical 
calculations of fireball dynamics that are computationally faster than hydrodynamic simulations. It is also a very 
flexible approach, that can be easily extended to include more complex situations, such as a continuous injection 
of energy at the reverse shock, and the sideways expansion in non-spherical ejecta. 

Some features of the numerical results for the remnant dynamical evolution are discussed and compared to the 
analytical results. We find that the ratio of the observer times when the jet half-angle reaches twice its initial value 
and when the light cone becomes wider than the jet is substantially smaller than predicted analytically. The effects 
arising from the remnant curvature and thickness further reduce this ratio, such that the afterglow light-curve break 
due to the remnant finite angular extent overlaps the weaker break due to the jet's sideways expansion. 

An analysis of the effects of the curvature and thickness of the remnant on the afterglow light-curves shows that 
these effects are important and should be taken into account for accurate calculations of the afterglow emission. 

Subject headings: gamma-rays: bursts - methods: analytical - radiation mechanisms: nonthermal 



1. INTRODUCTION 

Previous work on the dynamics of Gamma-Ray Burst (GRB) 
remnants was done analytically (Sari 1997, Vietri 1997, Wei 
& Lu 1998, Chiang & Dermer 1999, Rhoads 1999) or was 
based on hydrodynamical codes (Panaitescu & Meszaros 1998, 
Kobayashi, Piran & Sari 1999). Analytical treatments of the 
major afterglow features (Meszaros & Rees 1997, Vietri 1997, 
Waxman 1997, Wijers, Rees & Meszaros 1997, Meszaros , Rees 
& Wijers 1998, Sari, Piran & Narayan 1998, Rhoads 1999) have 
used asymptotic solutions where the remnant Lorentz factor is a 
power-law in radius and the post-shock energy density is deter- 
mined only by the ejecta bulk Lorentz factor. The advantages 
and drawbacks of the analytical and numerical approaches are 
evident: the former may not lead to sufficiently accurate results, 
but is more flexible and sometimes more powerful, while the 
second provides accuracy at the expense of substantial com- 
putational efforts, and may lack flexibility in including more 
sophisticated features. 

In this work we further develop the semi-analytical approach 
presented by Panaitescu, Meszaros & Rees (1998) (PMR98), 
with the aim of obtaining a method that allows fast computations 
of remnant dynamics, without making recourse to laborious hy- 
drodynamic numerical codes, while retaining enough physical 
details to allow an accurate calculation of the afterglow light- 
curves and spectra. 



2. REMNANT DYNAMICS 

2.1. Model Assumptions and Features 

We assume that, at any time during its evolution, the remnant 
is axially symmetric, i.e. there are no angular gradients in the 
shocked fluid, and all relevant physical parameters are functions 
only of the radial coordinate. This is equivalent to assuming that 
the initial energy distribution in the ejecta and the external gas 



are isotropic and that, in the case of a jet-like ejecta, the physi- 
cal parameters of the remnant respond on a short time-scale to 
the effect of sideways expansion of the shocked gas. We also 
assume that the Lorentz factor within the shocked fluid is con- 
stant. Other features taking into account in the model are: i) 
a delayed energy injection (Rees & Meszaros 1998), where for 
definiteness we shall use a power-law energy injection, ii) an in- 
homogeneous external fluid (Vietri 1997, Meszaros et al. 1998); 
for simplicity we consider an external density which is a power- 
law in the radial coordinate, which includes the homogeneous 
and the constant velocity pre-ejected wind cases, and Hi) the 
above mentioned sideways expansion, if the ejecta are not spher- 
ical (Rhoads 1999). Each of these features brings or modifies a 
specific term in the differential equations shown below, which 
give the evolution of the mass, kinetic and internal energy of the 
remnant. 



2.2. External Medium 

In the next sections we shall consider only homogeneous ex- 
ternal media, nevertheless the model we developed is applicable 
for a more general case, where the external gas density varies 
as a power-law with the radius: p ex (r) = Pd{r/rd)~ a , r& be- 
ing the deceleration radius, defined as the radius at which the 
fastest (initial) part of the ejecta, moving with Lorentz factor 
To, sweeps up an amount of external gas equal to a fraction F^ 1 
of its own mass Mq = E^/Fqc 2 : 



(3 - a)E 



1/3 



(1) 



Eq being the energy of the leading ejecta, i2o the solid angle of 
the cone within which it was released and rid the external particle 
density at r<j. The r& is the radius at which the deceleration of 
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the fireball due the interaction with the external gas becomes 
important. Note that has only a weak dependence on a. 

The continuous interaction with the external gas increases the 
remnant mass 



[dM] ex = D.(r) Pex (r)r 2 dr , 



(2) 



where £2(r) is the remnant solid angle at radius r, which is larger 
than the initial £2q due to the sideways expansion. The evolution 
of the half-angle 9 of the remnant is d6 = c s dt'/r, where c s = 
c I V3 is the comoving frame sound speed and dt' — (ficF) ~ 1 dr 
is the comoving time, F and (3 being the Lorentz factor and 
speed of the shocked ejecta. Using £2 = 2n(l — co&9), the 
remnant solid angle evolves as 



d£2 

dr 



/Q(4tt-Q) 

3(r 2 - 1) 



(3) 



The interaction of the remnant with the external medium de- 
celerates the remnant and heats it. Using energy and momen- 
tum conservation for the interaction between the remnant and 
the infinitesimal swept-up mass [dM] ea; , one obtains that the 
comoving frame internal energy of the newly shocked gas is 
(T — 1) times larger than its rest-mass energy, and the changes 
in the remnant internal energy U and total kinetic energy K are 
given by 

[dU] ex = A(F - \)c 2 [dM] ex , (4) 

[dK] ex = (Mc 2 + U) [dr] ex = - (F 2 - 1 )c 2 [dM] ex . (5) 

The multiplying factor A in the right-hand side of equation ( ij) 
was introduced to account for possible radiative losses in the 
shocked external gas. It represents the fraction of the internal 
energy of the shocked external gas that is not radiated away. 
Thus A = corresponds to a fully radiative remnant, which im- 
plies radiative electrons and strong coupling between electrons 
and protons, while A = 1 is for a fully adiabatic remnant. 

2.3. Delayed Energy Injection 

As suggested by Rees & Meszaros (1998) it is possible that 
the initial instantaneous deposition of energy in the fireball is 
not uniform in the entire ejecta, in the sense that some parts of 
it have been given more energy and have been accelerated to 
higher Lorentz factors, the slower ejecta catching up with the 
faster ones as these are decelerated by the interaction with the 
surrounding medium. The entire process is fully characterized 
by the energy distribution (dE/dF)i n j in the ejecta at the end of 
the initial phase of acceleration, with all other relevant quantities 
resulting from the kinematics and energetics of the "catching 
up". The velocity (3f of the incoming ejecta satisfies r — (3fct 
(free expansion), where t is the lab-frame time: t = td + 
/ [c/3(r)] _I dr, with (3(r) the speed of the decelerated ejecta. 
We assume that there is a significant delayed energy injection 
only at r > Td, so that equation ([!]) remains valid. Thus the 
Lorentz factor Tj of the delayed ejecta entering the decelerated 
part of the fireball is given by 



dr 



-f - 1 
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(6) 



Using the energy and momentum conservation, one can calcu- 
late the increase in the remnant kinetic and internal energies as 



result of the delayed injection: 

[dU] inj = [rT/(l - 00 f ) - 1] (?[dM] inj , (7) 

[dK} inj = (Mc 2 +U)[dF] mj = T / [l-r 2 ( 1-/3/3/)] c 2 [dM] mj , 

(8) 

where 



dE 

[dM] inj = ( — 



in j 



mL = F(T f ,WiNA, (9) 
Tf - 1 r 



is the infinitesimal injected mass, F(Tf,T) being a function 
that depends on the details of the delayed energy injection, 
and Mjnj the total mass of the delayed ejecta (see PMR98). 
In the numerical calculations we shall consider the particular 
case of a power-law injection: (dE /dT)i n j oc FJ S (Rees & 
Meszaros 1998) for r m < Tf < Td, where F m and Fd = 
r(r = Td) are the minimum and maximum Lorentz factors of 
the delayed ejecta. 

2.4. Adiabatic Cooling. Remnant Volume 

The delayed energy input at the reverse shock that moves 
into the incoming ejecta and the heating of the external fluid by 
the forward shock increase the internal energy of the remnant. 
This energy is lost adiabatically and radiatively. If they were 
acting alone, adiabatic losses would accelerate the remnant; in 
the presence of the external fluid, they re-convert internal energy 
into kinetic energy, which mitigates the remnant deceleration. 
This is described quantitatively by 

W\ad = -(7 - l)(drV'/V')U , (10) 



[dK] ad = (Mc 2 + U)[dT] ad = -T[dU] ad , (11) 

where V' — V^ s + Vp S is the comoving volume of the shocked 
ejecta (located behind the reverse shock) and of the swept-up 
external gas (behind the forward shock), and 7 is the adiabatic 
index of the remnant gas (7 = 4/3 for hot gas). In equation 
( |lO| ) d r V' denotes the infinitesimal variation of the comoving 
volume due only to the radial expansion of the gas, excluding 
the infinitesimal increases due to the addition of shocked fluid 
and to the sideways expansion (we neglect the adiabatic losses 
due to the sideways expansion, representing the acceleration of 
the outer parts of the fluid in the direction perpendicular to the 
radial direction of the flow). 

At this point we need a prescription for calculating d r In V . 
For the the radiative losses one also has to calculate the co- 
moving volume V' , to determine the comoving energy density, 
necessary for the computation of the magnetic field. We con- 
sider two models: Model 1, where we assume that the lab frame 
increase in the thickness of the shocked fluid (external or de- 
layed ejecta) is due only to the addition of new gas, and Model 2, 
where we assume that the comoving density of the two shocked 
fluids are uniform behind each shock, and have the values set by 
the shock jump equations. Only Model 1 is consistent with the 
assumption that F is constant in the entire shocked fluid, as this 
assumption implies that the remnant is neither dilating nor con- 
tracting in the radial direction. Model 2 implies the existence 
of a velocity gradient in the shocked fluid, nevertheless this gra- 
dient is expected to be small and the assumption of constant F 
within the remnant may still be safely used for the calculation 
of the relativistic effects. 
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2.4.1. Model 1 

The comoving volume is V = £2r 2 rA, where A is the lab- 
frame thickness of the remnant, determined by the relative mo- 
tion of the forward and reverse shocks: dA = ((3fs — /3fts)dt. 
The lab-frame reverse shock speed (3 R s can be calculated with 
a Lorentz transformation from j3' RS , the speed of the reverse 
shock measured in the frame of the ejecta entering this shock. 
The assumption that the frame thickness of the already shocked 
fluid remains constant (or equivalently that the flow velocity is 
uniform within the remnant) implies that d r In V' = dln(rr 2 ), 
therefore 

drV Mi„dr dT 



V 



= 2- 



and allows the calculation of j3' RS and (3fs from F, the remnant 
bulk Lorentz factor: 



RS 



/3T'[7(P- i) + r 



! (3fs 



(r-i)( 7 r + i) 
mr -i) + i] 



(13) 

where F' = FFf(l — (3(3 f) is the Lorentz factor of the shocked 
delayed ejecta measured in the frame of the incoming ejecta 
located just ahead of the reverse shock. After all the delayed 
ejecta has caught up with the remnant, or if there is no such 
delayed injection, (3rs is set equal to (3, the speed of the contact 
discontinuity. 



2.4.2. Model 2 

In this model, the two volumes V' RS and VLo can be calculated 
from the masses of the shocked gases (eqs. and [^]) and from 
the comoving densities, assumed to be uniform in each shell 
and having the same value as in the proximity of the shocks. 
The comoving densities p' RS and p' FS behind the reverse and 
forward shock are determined by the comoving densities of the 
un-shocked fluids, p'j and p ex , and by the Lorentz factors of the 
shocked gases, F' and F: 



Prs 



7r' + i 



7 . 



i 



p'f 



Pfs 



jF+1 



7 



1 



(14) 



Reverse Shock. The comoving density of the ejecta ahead 
of the reverse shock can be calculated by equating [dM]^ 
given by equation (^) with the mass r '/ pj£2/r 2 dZ swept up 
by the reverse shock as the remnant moves from r to r + dr, 
where dl = ct\d(3f\ — ((3f/(3 — l)dr is the infinitesimal 
lab frame distance relative to the contact discontinuity cov- 
ered by the reverse shock, and £lf (Ff, r) ^ £2 is the solid 
angle of the delayed ejecta. The assumption of isotropy 
within the ejecta beaming cone does not hold in the de- 
layed energy input case, our 1 -dimensional model including 
the delayed energy injection being applicable only to spher- 
ical ejecta for any lab-frame times and to beamed ejecta for 
times when the sideways expansion is negligible: £2/ ~ £2 ~ 
£2q. The end result is p'f = (F 2 e Pf/Q.or 2 ct)(dAI/dF) in j, so 
Vrs = (-^o + Jr [dM]i n j) / p' RS can be calculated by inte- 
grating equation ([)[) and using the first equation (|l4| ). From 
d r (\nV RS ) = — d(\np' RS ) and equations (Q) and (|14[), one can 
write 



d r V' 



RS 



ill 



Vrs 



dr 



ft 



l dr 



(15) 



where G(Ff,F) is a function of the details of the delayed 
injection. After the end of the delayed injection the comoving 
volume of the shocked delayed ejecta is considered constant. 

Forward Shock. The volume of the shocked external fluid 
V FS = (MoFq + C \dM] ex ) I p' FS can be calculated using 
the second equation (|l4| ) and the swept-up mass obtained by 
integrating equation (|^ From the second equation ( |l4| ) it can 
be shown that 



d r V FS Ml 

= 



dr 



7 



-rr + 1 



-dr 



(16) 



2.5. Differential Equations for Remnant Dynamics 

The remnant dynamics is given by the differential equations 
describing the evolution of the total kinetic and internal energies, 
coupled through the adiabatic losses: 

dV = [dU] inj + [dU] ad + [dU] ex , (17) 



dK = (Mc 2 + U)dF - [dK] inj + [dK] ad + [dK] ex , (18) 

where the quantities in the right-hand side terms are given 
by equations (§), (||, (@), (§), @i and (|l]). By substituting 
the term [df/] d from equation (|18[) in equation ( p7| ) ([dll] a d 
appears in the expression of [dK] a d - see eq. [|j"H]), one arrives 
at d[M (F- 1) + FU] =(Ff- l)[dM] inj , which simply states 
that the net variation of the total energy of the adiabatic remnant 
equals the input of energy through the delayed injection (global 
energy conservation). 

With the aid of all the relevant equations previously derived, 
equations ( fL7| ) and ( |T8| ) can be used to calculate dF/dr and 
dU/dr, i.e. the evolution of the flow Lorentz factor and of 
the co-moving internal energy. These equations are solved 
numerically together with equation (||) for Ff, equation #for 
£2, and the differential equation for the remnant mass resulting 
from equations (^) and (Q): 



dM s £2 

r— = F(F f ,F)M INJ + {3~a) 



dr 



, 3— a „ , T 

r \ M 



(19) 



The initial conditions for the set of differential equations for 
remnant dynamics are given by the values of the relevant quan- 
tities at r = r d (see PMR98): F f = F = 0.62r , U = 
0.62 AM Q c 2 , M = (1 + F-^Mq, £2 - £2 . 



3. NUMERICAL RESULTS FOR THE REMNANT DYNAMICS 

The remnant Lorentz factor F and the internal energy U de- 
termine the electron random Lorentz factor and the magnetic 
field, both necessary for the calculation of the afterglow emis- 
sion. Thus we are interested in solving the remnant differential 
equations to calculate the evolution of F and U with the observer 
time T, 



dT= (l+z)(l-f3)dt= (l + z) 



-ii*:, ( 2 o) 



where z is the source redshift. Equation ( |20| ) gives the time 
Tcd when the radiation emitted along the line of sight toward 
the observer and from the contact discontinuity arrives at Earth. 
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If most of the radiation comes from the fluid close to forward 
shock then it is necessary to calculate the observer time using 
the Lorentz factor of this shock. For a relativistic remnant 
rV s ~ V2T, thus 2W = T CD /2. 

Equations ((j]) and ( |l8| ) show that the remnant dynamics is 
determined by eq = E0/Q.0, the energy per solid angle in the 
ejecta, the jet initial solid angle, which determines when the jet 
sideways expansion becomes important, the parameters and 
a characterizing the surrounding medium, the remnant initial 
Lorentz factor To and the remnant radiative efficiency. In the 
case of an adiabatic remnant running into a homogeneous exter- 
nal medium, To cancels out from the expression for T(T), thus 
it is an irrelevant parameter, from the observer's point of view. 
However To is an important parameter for a radiative remnant, 
or if the external medium is not homogeneous. The remnant 
dynamics is also determined by the parameters of the delayed 
energy injection, which for a power-law injection are T m , To, s 
and Minj (or, equivalently, the entire injected energy Ejnj)- 
The effect of some of these parameters can be assessed from 
Figures 1-3. 

3.1. Spherical Remnants 

Figure 1 shows the evolution of — d log T/d log r for a spheri- 
cal remnant with no delayed energy input, running into a homo- 
geneous external medium. The non-relativistic phase, defined 
by r < 2, sets in at r < 10 rd for a fully radiative remnant and at 
r < 100 rd for an adiabatic one. An analytical treatment of the 
remnant dynamics leads to — d log T/d log r = (3 — a)/(l+A), 
as long as the remnant is relativistic. Thus, if a = 0, T oc r~ 3 
for a radiative remnant and T oc r~ 3 / 2 for an adiabatic one. 

1 /3 

These results hold for -C r <C T rd in the former case and 

2/3 

for rd <C r <C T ' ' in the latter. The values shown in Figure 
1 at early times (i.e. in the relativistic phase), are consistent with 
the analytical expectations. Due to the fact that the r~ 3 phase 
is short lived for a radiative remnant, this regime is not strictly 
reached for the case shown in Figure 1 (To = 500), where 
the steepest T-decay attained is oc r -2 . Only Lorentz factors 
To > 10 3 allow this phase to fully develop at very early observer 
times (T < 0. 1 day). In the case of a pre-ejected wind {a = 2), 
we obtained numerically the analytical results T oc r~'/ 2 and 
r oc r for an adiabatic and a radiative remnant, respectively 
(these cases are not shown in Figure 1). 

Figure 2 shows the effect of a delayed energy input on the 
dynamics of an adiabatic remnant, assuming a homogeneous 
external gas and a power-law distribution of energy per Lorentz 
factor in the delayed ejecta. The minimum Lorentz factor Y m of 
the ejecta determines the observer time when the injection ends. 
A sudden energy input (i.e. large parameter s), resembling the 
collision of a second shell with the leading fireball, may lead 
to a temporary flattening of T as a function of r, as shown by 
the small value of — d log T/d log r at T ~ 3 days for s = 10. 
The flux of the synchrotron radiation emitted by the remnant 
at a frequency v above the synchrotron peak v p (of vF v , the 
power-per-decade) is proportional to r 8+4/3 T 3 if the electrons 
radiating at v are adiabatic, and proportional to r 4+4/3 T 2 , if 
the same electrons are radiative, where (3 is the slope of the 
spectrum above v p : F v oc . For (3 ~ 1, as observed in most 
afterglows, the remnant flux varies like r 12 T 3 and r 8 T 2 for 
adiabatic or radiative electrons, respectively. This means that 
the afterglow corresponding to the remnant evolution shown in 
Figure 2 for s = 10 should exhibit a substantial brightening, 



with F v increasing as fast as T 3 (adiabatic electrons) or T 2 
(radiative electrons) at T ~ 3 days. 

3.2. Conical Remnants 

Before the effect of the sideways expansion becomes impor- 
tant, the bulk Lorentz factor of an adiabatic jet-like remnant 
running into and a homogeneous external medium is given by 

r = r rf (r/r d )- 3 / 2 . (21) 

The radius rj at which the remnant Lorentz factor has decreased 
to 6q , i.e. when an observer located on the jet's symmetry axis 
"sees" the jet's edge, if the jet sideways expansion until Tj is 
ignored, is given by 

rj = (T d e ) 2/3 r d . (22) 

Using equations ( |20| ) and ( |2l| ) it can be shown that we have 
r oc T 1 / 4 and T oc T~ 3 / 8 . For r > r, the area visible to the 
observer no longer increases as (rr) 2 oc T 5 / 4 , as it did when 
r < rj, but as (tOq) 2 oc T 1 / 2 , thus the light-curve decay should 
steepen by a factor T -3 / 4 . 

The remnant is still relativistic (r > 2) at rj if Qq ^ 30°. If 
the jet is sufficiently narrow, then the sideways escape may have 
an important effect on the remnant dynamics before the onset of 
the non-relativistic phase. Taking the radius at which the jet's 
half-angle is twice the initial one as the definition of the radius 
r b where the sideways expansion becomes important, and using 
the equation for rt, derived by Rhoads (1999), we obtain 

r b = (75/4) '/ 3 r 3 = [(75/8)r 2 2 ] 1/3 r d . (23) 

Equation ( ^3| ) is valid only if the remnant is still relativistic at 
r^,. Since the remnant Lorentz factor at r& is Ff, = (2/5^/3)6^ 
(Rhoads 1999), this condition reduces to 9 Q < 0.1 rad — 6° 
(£2 < 4 x 10~ 2 sr). 

Equations @ and @ show that r b /rj - (75/4) '/ 3 - 2.7, 

where we used F d — T0/V2, for consistency with Rhoads' 
results, which is close to the value 0.62Fo derived by PMR98. 
Therefore the jet edge effect should always be seen before that 
of the sideways expansion. Since Tar 4 for r < rt, the ratio 
of the observer times at which the sideways expansion and jet 
edge phases begin should be Tb /Tj — (n, / rj ) 4 ~ 50. 

As shown by Rhoads (1999), during the sideways escape 
phase r decreases exponentially with radius: 

r = r h e-( r -'' i >)/ r " . (24) 

The exponential constant can be cast in the form r e = 
(ToOo) 2 / 3 rd = 0.47 r b - Thus during the exponential regime 

- dlnr/d(r/r d ) = (T 9 )- 2 / 3 . (25) 

With the aid of equations ( p0| ) and ( j24| ) one can show that 
in the exponential regime T oc T -1 / 2 , thus the non-relativistic 
phase begins at T nr = (l/4)r 2 T b = (75 Ol)- l T b . Using equa- 
tion ([!]), the times Tj, TJ, and T nr can be calculated: 

1 /3 

507- = T b = (150 2 )T nr = 1.0 fi±M r [day] , 

(26) 
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where £0,54 is the initial energy per solid angle in units of 
1 54 erg sr" 1 , no is the external medium number density in cm -3 
andilo = 10 Qo,— 3 sr. 

The dynamics of adiabatic conical remnants is shown in Fig- 
ure 3, where we assumed a homogeneous external medium. 
As can be seen the exponential regime (i.e. the flattest part of 
each curve) is less evident for ejecta whose solid angle is larger 
than ~ 10~ 2 sr, when the onset of the non-relativistic regime 
occurs before the sideways expansion has a significant effect 
on the remnant dynamics. For a jet with Eq = 10 51 ergs and 
fto = 1(T 3 sr (0 O = 1°), which the case shown in Figure 3 
with dotted lines, equation (^) predicts that the jet edge is seen 
at Tj — 0.5 hours, the exponential regime starts at Tb = 1.0 
day and ends at T nr — 44 days. Numerically we obtain that 
T = 0" 1 atT_, = 1.0 hours, 8 being the jet half-angle, T b = 0.35 
days and T nr = 37 days for Model 1, and 1.2 hours, 0.35 days 
and 41 days, respectively, for Model 2. Note that the numerical 
and analytical results are in good agreement for T nr . 

The discrepancy for the Tj and Tb values arises from the fact 
that in the analytical derivation the effect of the sideways expan- 
sion on the remnant deceleration during the power-law phase 
was ignored. Because there is some sideways expansion during 
this phase, the jet half-angle 8 is increasingly larger than do, and 
r drops below 6~ x after it has reached the value Oq 1 , thus the 
analytical Tj underestimates the numerical one. Numerically 
we found that when r = 6~ 1 the jet angle is 8j ~ 1.2 9q. In the 
analytical treatment presented by Rhoads (1999) the increase in 
the swept-up mass due to the jet broadening during the power- 
law phase is ignored, which means that, for the same radius, the 
analytical F is larger than the numerical one, thus the analytical 
comoving time and the jet angle at given r are smaller than the 
values obtained numerically. Therefore the analytical Tb over- 
estimates the time when 8 — 2 6$. Numerically we found that 
at Tb given by equation ( |26| ) the jet angle is 8 ~ 2.5 do- 

For the jet whose dynamics is shown in Figure 3 (80 = 1°, 
To = 500), equation (£51) predicts that during the exponential 
phase — A\nT / d{r / r d) = 0.23, which is less than the values 
shown in Figure 3 at times after Tb and before T nr : 0.33 ± 0.08 
for Model 1 and 0.30 ± 0.07 for Model 2. This is consistent 
with the fact that — dlnr/d(?'/r,j) oc r~ l oc r^ 1 and that the 
numerical 7% is smaller than the analytical one. 

The same conclusions have been reached by comparing the 
analytical and numerical results for other values of £Iq and To, 
the most important them being that Tb/Tj = 7 -j- 10, which is 
5-^7 times smaller than obtained analytically. 

4. SYNCHROTRON LIGHT-CURVES FROM BEAMED EJECTA 

The calculation of the afterglow light-curve is described in 
PMR98, and consists in integrating the remnant synchrotron 
emission over its dynamical evolution, over the electron distri- 
bution in each infinitesimal layers of swept-up external gas and 
over the angle relative to the jet axis. The electron distribution is 
initialized as a power-law and evolved through synchrotron and 
adiabatic cooling. We assume an electron index p = 3, and that 
electrons acquire 10% of the internal energy after shock accel- 
eration. We also assume a turbulent magnetic field which stores 
10 -4 of the shocked gas internal energy, i.e. the magnetic field 
intensity is two orders of magnitude weaker than the equipar- 
tition value, which leads to negligible radiative losses and an 
adiabatic remnant evolution. For a more general case where 
there could be substantial radiative losses, some of the remnant 
dynamics equations are altered as following. The factor A in 



equation (|4j) is set to 1 and a new term is added to equation (|17|), 
representing the radiative losses. The second term in the right 
hand side of equation (|l7|), representing the adiabatic losses, 
is split into a term for the adiabatic losses of the baryons and 
another one for the electrons. The electronic adiabatic losses 
act simultaneously with the radiative losses and are put together 
in a differential equation for electron cooling (equation [13] in 
PMR98). 

The calculation of the photon arrival time takes into account 
both the geometrical curvature of the remnant (photons emitted 
by the fluid moving off the observer's line of sight toward the 
center of the shell arrive later than those emitted along this cen- 
tral line of sight) and the thickness of the shocked fluid, whose 
light crossing-time is not negligible compared to the observer 
time since the main burst. Here we assume that there is a neg- 
ligible mixing behind the forward shock, such that the electron 
distribution in any remnant infinitesimal layer is the one injected 
when that layer was added to the shocked structure, evolved 
through adiabatic and radiative losses. The shell thickness at 
the time when a new infinitesimal layer of swept up external gas 
is added can be obtained from the remnant comoving volume, 
whose calculation is described in §2.4. 

The effect of the sideways expansion on the optical afterglow 
seen by an observer located on the jet axis is shown in Figure 4a. 
The initial energy per solid angle is the same for all remnants, 
only the jet initial half-angle do is changed. The afterglow 
brightness should be independent of 9q until Tj, when the flow 
Lorentz factor has become sufficiently low that the observer sees 
the edge of the jet. This feature is better seen if the sideways 
expansion is "switched off", because in the case where it is taken 
into account there is a non-negligible jet broadening until Tj. 

For the afterglows shown in Figure 4a F decreases to 6~ l 
at Tj = 0.005, 0.086, and 1.8 days for 6 = 1°, 3°, and 9°, 
respectively, if the jet broadening is not taken into account 
(i.e. 8 = 8q st all times), and 0.009, 0.16, and 4.1 days, re- 
spectively, if the sideways expansion is accounted for. These 
are the times when photons emitted from the forward shock 
along the remnant center-observer line (which is the jet axis in 
this case) arrive at the observer. Photons emitted from the fluid 
located closer to the contact discontinuity arrive up to twice 
later. For T < Tj, photons emitted from the forward shock 
regions moving at an angle r _1 off this central line of sight 
arrive at T = (1 — cosF~ l )(r /c) ~ (2F 2 )~ l (r /c), which is 
factor 8 larger than the arrival time from the forward shock 
Tfs = ( 16r 2 ) ~ 1 (r/c), as can be shown using equations ( ^o| ) 
and (|l]). 

The times when the half-angle of the jets whose afterglows are 
shown in Figure 4a reach twice their initial values are = 0.08, 
1.3, and 27 days, for 8 = 1°, 3°, and 9°, respectively. The optical 
light-curves shown in Figure 4a steepen smoothly around X&, 
while the light-curves of the non-broadening jets maintain the 
decay slopes they had before Tb (of course, Tb has no meaning 
for a jet of constant opening). The light-curve steepening that 
can be seen for the non-broadening 80 = 9° jet around T = 10 
days is due to the passage of the cooling break through the 
optical band. 

It can be noticed that the slopes of the light-curves for non- 
broadening jets shown in Figure 4a are not constant after the 
Tj's given above and before Tb, as the remnant geometrical 
curvature delays the photon arrival time from regions off the 
jet axis. Moreover, the received power per solid angle being 
proportional to [r(l — (3 cos <5)] -4 , where S is the angle relative 
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to the central line of sight at which an infinitesimal emitting 
region moves, implies that at Tj this power per solid angle from 
the jet edge (S = 6j = FJ 1 ) is 16 times smaller than that from 
(5 = 0, which leads to the conclusion that it should take longer 
than just Tj to see an afterglow dimming rate in excess of that 
existent until ~ Tj (this is confirmed by the numerical results, as 
shown below). Put together with the fact that T b ~ (7 10)1} 
(see previous section) this suggests that the effects arising from 
seeing the jet edge and from the sideways expansion may not 
be so clearly distinguishable for the observer. 

For the sideways expanding jets shown in Figure 4a, the non- 
relativistic phase begins at T nr = 11, 18 and 32 days (forward 
shock times) for 9q = 1°, 3°, and 9°, respectively, while for the 
jets where sideways expansion was not taken into account in the 
dynamics T nr = 59 days, independent of 8q. 

The effect of the sideways expansion in the case where the 
observer is located off the jet axis, at an angle 8 obs relative to 
this axis, is illustrated in Figure 46. The major difference from 
the 6 obs = case is that, shortly after the light-curve rises, 
which happens later for larger 6 obs , the broadening of the jet 
yields a brighter afterglow than in the case where the sideways 
expansion is ignored. This is due to the fact that for broadening 
jets there is some shocked fluid approaching the observer line 
of sight toward the remnant center, along which the relativistic 
effects are maximal, while for non-broadening jets it is only the 
decrease of the remnant bulk Lorentz factor that "brings" the 
observer into the cone of the relativistically beamed radiation 
and thus to see the afterglow. 

In Figure 5 we illustrate the importance of integrating the 
remnant emission over the angle relative to the observer and 
by taking into account the finite thickness of the source, i.e. the 
time delays introduced by the location within the remnant where 
the radiation is released. The curves shown with thick broken 
lines correspond to afterglows where both of the above effects 
are taken into account. The curves shown with thin broken 
lines represent afterglows where the integration over the polar 
angle was "switched off", in the sense that the light-curve has 
been calculated as if all the emitting fluid were moving directly 
toward the observer. However, to account for the fact that at 
times T < Tj and in the presence of the relativistic effects the 
observer receives radiation only from the fluid moving within 
r _1 off the central line of sight, the afterglow's brightness at 
T < Tj has been corrected by a factor (r6>)~ 2 ~ (r6> )~ 2 . As 
can be seen, ignoring the effect arising from the geometrical 
curvature of the source leads to a substantial overestimation of 
the early radio afterglow (Figure 5a), introduces errors of more 
than 1 mag in the early optical afterglow and of ~ 1/2 mag at 
later times (Figure 56), and underestimates the X-ray afterglow 
by a factor < 4. 

The evolutions of the slope of the afterglows obtained by ac- 
counting for the shell curvature are given in Figure 5d, where 
we also show with vertical lines the important forward shock 
times. For the parameters chosen here for the electron frac- 
tional energy and for the magnetic field intensity, the spectral 
cooling break is always between optical and X-ray and above 
the synchrotron peak, thus the electrons radiating at this peak, in 
radio and in optical are adiabatic, while those radiating in X-ray 
are radiative. Following the analytical treatment of Meszaros et 
al. (1998), one can derive the following asymptotic afterglow 



light-curves for the case discussed here: 
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where p — 3 for the light-curves shown in Figure 5. The 
afterglow slope is expected to steepen at all frequencies by —3/4 
when the jet edge is seen, while the changes in the same slope 
due to the jet sideways expansion are expected to be — 1/ 12 in 
radio, — p/4 = —3/4 in optical, and — (p — l)/4 = —1/2 in 
X-ray. It should be noted here that in the derivation of equations 
(p7[)-(p9|) the time evolution of the remnant specific intensity 
at the synchrotron peak was calculated assuming that all the 
electrons radiate at this peak (the electrons being adiabatic), 
which is a rough approximation given that the injected electron 
distribution (and implicitly the peak energy of their synchrotron 
emission) changes substantially during an adiabatic timescale, 
which is comparable to the age of the remnant. 

The afterglow slopes shown in Figure 5d for the case where 
the "angular effect" is ignored exhibit the trends predicted by 
equations (p7|)-(p9|). The largest difference is shown by the 
optical light-curve at times between ~ Tj and until after T b . 
It results mostly from the fact that at these times the optical 
band is above but close to the synchrotron peak, the integration 
over the electron distribution in the shell leading to a spectrum 
that does not exhibit a pure — (p — 1 ) /2 slope (the electrons 
radiating in optical are adiabatic), as assumed in the derivation 
of the analytical light-curves, but to a smooth transition between 
the 1 /3 slope below the synchrotron peak and the — p/2 slope 
above the cooling break (see Figure 6). 

Note from Figure 5d that taking into account the curvature of 
the emitting shell delays the jet edge break, as expected, until 
times comparable to T b . The effect of the remnant thickness 
on the afterglow slope evolution (not shown in Figure 5d) is 
similar, but weaker, to that of the remnant curvature. Figure 6 
shows that the spectrum calculated by taking into account the 
time delays and Lorentz boosting factors over a curved shell is 
harder than if the shell is assumed planar. A photon emitted 
by the gas moving off the observer's central line of sight is 
at most twice less blueshifted than one emitted along this line, 
nevertheless the former was radiated at a smaller radius, when F 
and the synchrotron peak frequency (which scales as T 4 ) where 
larger, the end result of the integration over the equal arrival 
time surface being a harder spectrum. 

The upper three panels of Figure 5 also show the light-curve 
obtained with Model 1 if it is assumed that all the radiation 
is emitted from the location of the forward shock, i.e. if the 
light travel time across the source is neglected. The differences 
are minor in the X-ray because the electrons emitting at high 
energies are radiative and cool fast, thus the X-ray emission 
from the fluid that is not located close to the forward shock is 
negligible, but significant differences can be seen in the radio 
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and optical light-curves, as the electron radiating in these bands 
are adiabatic, and thus occupy the entire volume of the remnant. 
Also shown in Figures 5a and 5b with continuous lines are the 
afterglows calculated with Model 2 for the remnant comoving 
volume. Note that the largest difference between the optical 
light-curves within two models amounts to ~ 1 mag in the late 
afterglow. These differences are due mostly to the prescription 
for the comoving volume, which leads to different magnetic 
fields and shell thicknesses, and not to the adiabatic losses. 

5. DISCUSSION 

We have presented an analytical treatment for the dynamics of 
an expanding fireball, capable of following its evolution from 
the onset of the deceleration phase (r ~ r<j) until arbitrarily 
large times. The differential equations for the remnant dynam- 
ics given here are valid in any relativistic regime. The major 
assumption underlying the analytical derivations is that, at any 
time, the remnant is axially symmetric. 

This analytical treatment takes into account a possible de- 
layed energy input resulting from an impulsive but uneven de- 
position of energy in the ejected material. For beamed ejecta, 
it also takes into account the intensification of the remnant de- 
celeration due to the increase of the solid angle of the remnant 
and, thus, of the rate at which it sweeps up external gas. The 
results presented in the previous section illustrate the effect of 
these two factors. The treatment of the adiabatic losses and 
the calculations of the magnetic field and remnant thickness re- 
quire a prescription for how to calculate the remnant volume. 
We considered two models for this: Model 1 is based on the as- 
sumption that, if the accumulation of swept-up gas is subtracted, 
the remaining increase of the lab frame volume is due only to 
the r 2 increase of the remnant area, and Model 2, which is based 
on the assumption that the density profile behind each shock is 
uniform. The two models for the comoving volume calculation 
lead to significant differences in remnant dynamics when there 
is a sharp delayed energy input, as shown in Figure 2, and in 
the case of beamed ejecta, as shown in Figure 3. There are also 
significant differences among the radio and optical afterglows 
calculated within the two models for the remnant volume, as 
illustrated in Figures 5a and 5b. 

We have compared the features of the dynamics of beamed 



ejecta calculated numerically with some analytical results de- 
rived from those presented by Rhoads (1999). The most impor- 
tant difference is that the ratio of the observer times T),, when 
the jet half-angle 6 reaches twice its initial value due to the side- 
ways expansion, and Tj, when the remnant Lorentz factor has 
decreased to 9 , is analytically overestimated to ~ 50, while 
the numerics lead to a value ^10. Tj is the time when the after- 
glow is expected to steepen by T~ 3 / 4 due to the finite angular 
opening of the ejecta, while at Tj, the afterglow should further 
steepen by at most T~ p / 4 , where p is index of the injected 
electron distribution, due to the intensification of the remnant 
deceleration caused by the continuous jet broadening. The af- 
terglow slope evolutions shown in Figure 5d suggest that, due 
to the remnant curvature, the Tj defined above actually under- 
estimates the time when the finite angular extent of the remnant 
yields a substantial steepening of the light-curve, and that jet 
edge break overlaps with the weaker light-curve steepening due 
to the jet sideways expansion. 

We have compared the numerical light-curves with the 
asymptotic behaviors expected analytically and we have ana- 
lyzed the importance of integrating the remnant emission over 
its angular opening, as the Doppler boosting and the photon 
arrival time are functions of the angle between the direction 
of motion of the emitting fluid region and the direction toward 
the observer. Figure 5 shows that this effect is quite important, 
and should be taken into account for accurate afterglow calcu- 
lations. We have also found that the thickness of the source is 
an important factor, though less than the shell curvature, in the 
calculation of radio and optical afterglows. Lastly, the after- 
glow spectrum at frequencies around the peak and the cooling 
breaks, and between these breaks (if they are not too far from 
each other), may be poorly approximated as a power-law, par- 
ticularly if the electrons radiating at the frequency of interest are 
adiabatic (see Figure 6), implying that they have been injected 
over a wide range of radii and bulk Lorentz factors. In this case, 
a numerical integration of the electron injection and cooling is 
necessary to obtain an accurate afterglow emission. 
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Fig. 1. — Evolution of — dlogT/dlogr for adiabatic (A = 1) and radiative (^4 = 0) spherical remnants, with no delayed injection 
and homogeneous external medium. Parameters: z = 1, To = 500, rid = 1 cm~ 3 and eo = E0/Q.0 as given in the legend, in units 
of erg sr _1 . Thick curves correspond to Model 1, while thin lines are for Model 2. Obviously, the numerical results are the same for 
both models of adiabatic losses if the remnant is fully radiative ( A = 1). Note that a larger energy per solid angle in the ejecta leads 
to a longer relativistic phase. 




Fig. 2. — Effect of a power-law delayed energy input on the dynamics of an adiabatic spherical remnant interacting with a homogeneous 
external medium. The inset shows the evolution of the kinetic and lab frame internal energy, in units of initial remnant energy. The 
parameters of the injection are F m = 10, Einj — Eo and s is given in the legend. Other parameters are: z = 1, £0 = 10 52 erg sr~', 
rid = 1 cm~ 3 . As for Figure 1, Model 1 solutions are shown with thick lines, while Model 2 solutions are indicated with thin curves. 
The continuous thin curve is for the case with no delayed injection (added for comparison). Note that the differences between the 
solutions obtained with the two models for adiabatic losses are larger at times when most of the injection takes place, indicating that 
the adiabatic cooling of the delayed shocked ejecta is the source of these differences. 
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Fig. 3. — Evolution of — d In F/ dr for adiabatic beamed ejecta. The remnant dynamics depends on the initial solid angle of the ejecta 
(given in the legend in steradians) and on the initial energy per solid angle (given in the legend in erg sr _1 ). Other parameters: z = 1, 
rid = 1 cm" 3 , a — Q,T() = 500. The meaning of thick and thin curves is the same as for Figures 1 and 2. The parts indicated with a 
continuous curve for the (£lo = 10~ 3 , eo = 10 54 ) models correspond to T between T b and T nr , i.e. when the remnant opening angle 
is more than twice the initial one and the remnant is still relativistic (r > 2). 
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Fig. 4. — Optical light-curves and the effect of sideways expansion for beamed ejecta. For both panels the R c magnitudes are 
calculated with Model 1. Thick curves are for the light-curves obtained by taking into account the sideways expansion of the ejecta, 
both in the dynamics and the integration of light, while thin curves are for the case where the jet broadening is ignored. Parameters: 
z = 1, rid = 1 cm~ 3 , a = and eo = 10 53 ergsr _1 . The electron and magnetic field energy densities are 10 _1 and 10~ 4 of the 
internal energy density. The distribution of the shock-accelerated electrons is assumed to be a power-law of index —3. Panel a: the 
observer is located on the axis of the jet, whose initial half-angle is given in the legend. Note that the jet progressive broadening leads 
to dimmer afterglows. Panel b: a jet of initial half-angle 9 Q = 3° is seen at different angles. Note that for 9 obs > 9 Q the sideways 
expansion yields afterglows that have brighter peaks than in the case where the jet broadening is not taken into account. 
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Fig. 5. — Effects of the remnant geometrical curvature and thickness on the observed light-curve. The upper three panels show the 
Model 1 radio, optical and X-ray afterglows for a jet with &o = 1° seen by an observer located on its symmetry axis. Other parameters 
are as for Figure 4. The light-curves whose calculations took into account both the source curvature and thickness are shown with 
thick broken lines, those for which the curvature was ignored are represented with thin broken lines, while stars are for the case 
where it was assumed that all the radiation is emitted from the location of the forward shock. The evolution of the radio (dotted line), 
optical (dashes), and X-ray (long dashes) light-curve slopes for the first two cases are shown in panel d. The relevant forward shock 
times are indicated on the abscissa. To = 0.005 days and Tr = 3.1 days are the times when the synchrotron peak passes through the 
optical and radio bands, respectively. For comparison, panels a and b also show with continuous thin lines the light-curves obtained 
with Model 2. The X-ray afterglows in the two models for the comoving volume are indistinguishable. 
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Fig. 6. — Spectral evolution of the afterglow shown in Figure 5 and the effects of the remnant curvature on the observed spectrum. 
Thick dotted curves are for spectra whose calculation included the shell curvature, while thin dotted lines are for the case where it was 
ignored. The squares connected by a continuous line indicate the synchrotron peak v p of the radiation integrated over the power-law 
distribution of the electrons freshly injected by the forward shock, while circles show the peak v c of the radiation from the electrons 
that cool radiatively on a timescale equal to that of the adiabatic losses. In both cases the Doppler factor corresponding to the motion 
toward the observer (the jet axis in this case) was used. From right to left the observer times for each spectrum are: 0.0 1 , 0. 1 , 1 , 10, 
and 100 days. Note that the integration over the equal arrival time surface yields harder spectra and spectral peaks above v p . Also 
note the smooth transition from the 1 /3 power-law spectrum at v < v v to the — p/2 = —3/2 one at v > v c , and the absence of a 
perfect — (p — l)/2 = —1 power-law for v p < v < v c . 



